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Abstract. We present results about the distribution and 
kinematics of the molecular environment of the AGN in 
NGC 1068, over a 1.5" x 3.5" region around the central 
engine in NGC 1068, derived from H2 line emission de- 
tected with ISAAC at VLT/ANTU on ESO/Paranal. The 
H2 emitting molecular gas is found to be distributed along 
the East- West direction and with two main peak emission 
(knots) located at a distance of about 70 pc from the cen- 
tral engine. The eastern H2 knot is more intense than the 
western one. The line profiles mapped across the entire 
1.5" X 3.5" region, at a spatial resolution of 0.3" x 0.45", 
appear to be quite complex with either a blue or red wing. 
At first order, we find a velocity difference of 140 km/s 
between the two knots; if interpreted as quasi-keplerian 
velocity, this implies a central enclosed mass of 10* M©. 
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1. Introduction 

Given the proximity of NGC 1068 (14.4 Mpc and corre- 
sponding scale of 70 pc per 1") and the predicted size of 
the molecular torus in an Active Galactic Nucleus (AGN) 
-from 1 to 100 pc-, any trace or signature of a molecular 
torus in NGC 1068 must be searched from data collected 
under sub-arcsec image quality. Several recent discoveries 
point towards the presence of a conspicuous and struc- 
tured molecular/dusty environment around the central en- 
gine of NGC1068 (e.g. Gallimore et al 1997, Rouan et al 
1998, Marco & AUoin 2000, Schinnerer et al 2000, and ref- 
erences given in these papers). One piece of information 
still missing to ascertain the existence of a rotating torus 
-in addition to other possible molecular components-. 
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is the kinematical status of the molecular/dusty mate- 
rial. Such information is available for the cold molecular 
gas, from recent interferometric work in the CO (2-1) line 
(Schinnerer et al 2000): the authors infer the presence of 
a warped disc of cold molecular gas. In this study we have 
chosen to probe the molecular gas through the II2 1-0 S(l) 
line (rest A ~ 2.12/im). Being related to the hot molecu- 
lar gas and having specific excitation mechanisms, it pro- 
vides complementary information to that derived from the 
CO transition which traces the cold molecular gas. After 
the pioneering discovery of the II2 2.12 /im line emission 
in NGC 1068 by Thompson et al (1978), a first attempt 
to image the AGN of NGC 1068 in the H2 2.12 fj.m line 
has been reported by Rotaciuc et al (1991), covering a 
10" X 10" region at a resolution of « 1". It is imperative 
to push further the spatial resolution. 

Adaptive optics high resolution K, L and M band im- 
ages of the AGN have unveiled the presence and structure 
of hot to warm dust (Rouan et al 1998, Marco & AUoin 
2000) within the 1" x 1" region around the central engine. 
We have selected from these observations two particular 
directions for our kinematical study: (a) PA = 102° pos- 
sibly tracing the equatorial plane of the molecular /dusty 
torus around the central engine; (b) PA = 12° which fea- 
tures the axis of the torus and is found to be close to the 
axis of the ionizing cone to the North-East. The location 
of the central engine - only visible directly in the IR as 
an unresolved core carrying around 90% of the emission 
and in the radio as the radio source SI in Muxlow et al 
(1996) - is taken here as that derived by Marco et al (1997) 
from simultaneous K and I band high resolution imaging 
with adaptive optics. Given the error bar on this position 
(±0.15") , it is coincident with the astrometric position 
of the 12.4 fj,m unresolved core (Braatz et al 1993), the K 
band peak observed by Thatte et al (1997) and the center 
of symmetry of the polarization pattern in the near-IR and 
mid-IR (Lumsden et al 1999). We leave aside attempts at 
locating the central engine from UV data which are quite 
sensitive to dust extinction and provide so far discrepant 
results (Capetti et al 1995a & b, Kishimoto 1999). 
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We present and discuss results of long-slit spectroscopy 
in the near-IR, obtained with ISAAC at VLT/ANTU on 
ESO/Paranal. In this Letter, we concentrate on the H2 
2.12 /im line emission observed over a 1.5" x 3.5" region 
through the central engine at PA = 102°. The kinematics 
of the hot molecular gas in particular is investigated to 
test the presence of an eventual super massive object in 
NGC 1068. 

2. Data collection and Reduction 

The observations were performed using the SWSl short 
wavelength arm of the instrument ISAAC attached to the 
Nasmyth focus of ANTU (Moorwood et al 1999, Cuby et al 
2000). The measured seeing value was of 0.5" (FWHM) in 
the K band. The slit width was set to 0.3" while its length 
was of 2' at PA = 102°. A spectrum was first obtained with 
the slit centered on the near-IR unresolved core, imaged 
prior to the spectroscopic observations. Then the position 
of the slit was offset by 0.3" and 0.6" on each side of the 
core, to the North and to the South, providing a complete 
mapping of a 1.5" x 1' area. A complete description of the 
data collection and reduction procedures is given in Gal- 
liano & AUoin (2001). The final spectral resolution at the 
observed 2.15 /im wavelength of the H2 line is 35 km/s. 
From the reduced 2D spectra, we extracted a series of ID 
spectra, through a window 3 pixels-high (i.e. 0.45") along 
the slit, and with a sliding step of 1 pixel (0.15"). In the 
very central area, the intense continuum produces a fring- 
ing pattern at the level of 4% (peak to peak) that cannot 
be fully corrected, leaving some residuals which prevent 
the measurement of eventual faint II2 line emission at a 
distance less than 15 pc from the central engine. 

3. Results and Discussion 

3.1. H2 line profiles 

We provide in Figure ||, a set of the H2 line profiles ob- 
served at various positions across the 1.5" x 3.5" central 
area extended along the equatorial plane of the suspected 
torus (PA = 102°). At each position, the displayed line 
profile corresponds to an 0.3" x 0.45" emitting patch. We 
did not deconvolve the individual successive spectra: given 
the seeing (0.5" in the K band), each spectrum is only 
moderately contaminated by its neighbours. The H2 line 
profiles are presented in the following manner: for each 
frame, offset coordinates of the patch to which the line- 
profile refers are given with respect to the central engine 
location. The horizontal axis provides the offsets along 
PA= 102° and the vertical axis the offsets along PA= 
12°. In each frame, the black (thin) line features the line 
profile observed to the West of the AGN (from —1.5" to 
—0.3"), while the red (thick) line features the line profile 
observed to the East (from +0.3" to -1-1.5"). Therefore, 
each frame provides the line profile of a given patch and 
the line profile of its symmetrical counter part with respect 



to the axis of the suspected torus. Because of a large dif- 
ference between the H2 line fluxes to the East and West of 
the central engine, we have rescaled by a factor 3, in Fig- 
ure]^, the line profiles detected to the West. In most of the 
frames, the H2 line shows a complex profile: from double 
peak (North) to the presence of an extended wing ( blue 
wing to the West and red wing to the East). The veloc- 
ity difference between the main peaks in two symmetrical 
patches is maximum for patches [0";0.9"] and [0";— 0.9"] 
and reaches a value of AV = 140 km/s. 

3.2. Spatial distribution of the H2 line emission 

For each 0.3" x 0.45" emitting patch we have measured 
the total integrated H2 line emission and, from these val- 
ues, we have reconstructed the 2D map in H2 line emis- 
sion shown in Figure ^. The existence of two symmetrical 
emission knots (eastern and western knots) with large flux 
difference is conspicuous. This map can be compared to 
the early image in the II2 line provided by Rotaciuc et 
al (1991, their figure 2), within the limitation of differ- 
ent spatial resolutions and spatial coverage. Indeed, the 
ISAAC reconstructed H2 map is consistent with the in- 
ner part of the image by Rotaciuc, although the more 
limited area covered by the ISAAC data set may miss 
the maximum of the western knot (which is also slightly 
shifted to the South). A detailed flux comparison between 
H2 data sets obtained previously (references in Section 1) 
and the ISAAC data set is presented in Galliano & AUoin 
(2001). We concentrate here on the positional/kinematical 
aspects measured for the first time from this ISAAC II2 
data. The strong eastern H2 knot, located about 1" to 
the East of the central engine, is also well identified in 
velocity space (Figure |l|). In addition to the eastern and 
western knots, extended and asymmetric wings on the line 
profiles suggest the presence of an extended source of II2 
emission, possibly in the form of an outflow (Galliano & 
AUoin 2001). We have also reconstructed a 2D map of the 
average continuum between 2.1^m and 2.2/xm (not shown 
here) from which we have positioned the cross featuring 
the central engine in Figure 2. 

3.3. Kinematics of the H2 molecular material 

As a first order result, a jump in radial velocity of about 
140 (±5) km/s is detected between the eastern and west- 
ern knots at 70 pc from the central engine. The [Position- 
Velocity] diagram (Figure ^ cnlights this result. If the 
material associated with these knots pertains to a struc- 
tured ensemble of molecular clouds in keplerian orbit, this 
implies the presence of an enclosed object with dynamical 
mass 10* Mq. A comparable value is found by Schinnerer 
et al (2000) from CO line emission observed on a scale of 
±70pc. Yet, the enclosed central mass derived from H2O 
maser observations on a scale of less than 1 pc, is only 
1.5 X 10^ Mq (Greenhill & Gwinn 1997). It should be no- 
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Fig. 1. Position- Velocity diagram for the spectrum at 

PA=102° through the central engine. The continuum was 
fitted and substracted. The contours are at 10%, 20%,..., 
100% of the eastern peak. Inside the central marked box, 
the residual fringing pattern has been removed for clarity 
of the diagram. Outside the box, the fringing is absent. 

ticed as well that the rotation axis of the molecular and 

mascr "pseudo" discs differ by at least 30°. This result 
suggests that the mass enclosed within the 70 pc radius 
region around the central engine in NGC 1068 is shared 
among different components. 

At second order, across the western knot itself, an in- 
crease of the H2 line velocity is observed towards the cen- 
ter, as expected in the case of pseudo-kcplcrian rotation of 
the system, but a similar behaviour does not appear in an 
obvious way in the eastern knot where the brightest peak 
seems to remain at constant bhicshift when approaching 
the center. Considering the large flux difference between 
both sides, we cannot exclude a superposition of compo- 
nents along the line of sight to the East of the central en- 
gine, or the existence of peculiar kinematics as discussed 
from CO line data by Schinnerer et al (2000). A more 
refined analysis of the velocity field of the H2 emitting 
molecular material is in progress, aiming in particular at 
fitting the line profiles across the entire region. 



is detected as well in CO, at the same location. On the 
western side, the ISAAC spatial coverage does not allow 
to check whether the CO peak - located 1.5" to the West 
and 1" to the South - is also emitting in the H2 line. 
However, we do detect in H2 the tip of the western CO 
knot. The spatial coincidence of CO and H2 emission in 
the eastern and western knots is interesting and indica- 
tive that, in these knots, H2 is emitted in photo (or/ and 
X-ray) dissociation regions (FDR or XDR) at the surface 
of molecular clouds. Given the AGN environment, rich in 
UV photons, X-rays and shocks, it is still unclear what 
would be the main source of excitation of the H2 line. For 
the time-being, let us simply assume that the H2 line emis- 
sion originates from molecular gas at temperature around 
2000 K (e.g. Hawarden et al 1995): we find that the total 
observed H2 line emission in the 3.5" x 1.5" central area 
is of (20±6).10~^^ erg.s~^.cm~^, corresponding to a total 
luminosity of 1.3 x 10^ Lq. Following Thompson (1978), 
an amount of 3500 Mq of hot molecular gas is needed 
to produce this emission. Conversely, CO line emission is 
thought to arise from cold molecular gas with tempera- 
ture at least one order of magnitude less (e.g. Barvainis et 
al, 1997). The total mass of molecular hydrogen derived 
from the CO line flux by Schinnerer & al (2000), over the 
same region, is of the order of 10'' Mq, 4 orders of magni- 
tude larger than the hot emitting H2 molecular gas. Such 
a result is not surprising and is consistent with the PDRs 
or XDRs origin of H2 . Still to be understood is the more 
diffuse H2 emission which is traced in the wings in the II2 
line profiles (Figure 1) and we are developping a kincmati- 
cal model for this extended component (Galliano & AUoin 
2001). 

Another interesting point to resolve is whether the dy- 
namical mass 10^ M0 enclosed within a 70 pc radius re- 
gion around the central engine (as derived both from cold 
and hot molecular gas kinematics independently) corre- 
sponds to a super massive confined object or a mass- 
distributed component (stellar cluster and molecular /dust 
gas). As there is most probably a contribution from each 
of these components, their respective share is an issue we 
intend to address in the near future (Galliano & AUoin 
2001). 



3.4- Com,parison with CO high resolution map and 

Discussion 

The integrated H2 line map derived from ISAAC data can 
be compared with the ^^C0(2 — 1) map recently obtained 
at TRAM with an 0.7" x 0.7" beam by Schinnerer et al 
(2000). The strong H2 line emission of the eastern knot 
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Fig. 2. The H2 line-profile data: 
each frame corresponds to an 0.3" x 
0.45" emitting area. The center of 
the area is offset from the central en- 
gine position [0" x 0"] by the quanti- 
ties provided along the right side of 
the Y-axis (offset along PA = 102 °) 
and the upper X-axis (offset along 
PA = 12°), the X-axis offset being 
positive to the East and negative to 
the West, virhile the Y-axis offset is 
positive to the North and negative 
to the South. The lower-left frame 
provides as well the velocity scale 
applicable to all frames. The thin 
(black) line refers to the West of 
the AGN, while the thick (red) line 
refers to the East. As the first aim 
is to show the velocity shifts, the 
western profiles have been rescaled 
in flux by a factor 3. The flux scale 
is given on the left side of the Y- 
axis: the upper (red) one refers to 
the eastern profiles, the lower (black 
one) to the western ones. 
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Fig. 3. Restored map of the 1.5" x 
3.5" central area, in the H2 1.12 
fim integrated line emission. The in- 
tensity scale is linear and reaches 
2 X 10~^''erg.s~^.cm~^. The cross lo- 
cates the maximum of the contin- 
uum emission as found from this 
ISAAC data set and features the 
location of the central engine as 
derived from K band adaptive op- 
tics measurements by Marco et al 
(1997). 
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